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Abstract 

We study model-independently the implications of non-standard scalar and pseudoscalar in- 
teractions for the decays b — ► sj, b — > sg, b — ► s£ + £~ (£ = e, fj,) and B s — > fi + fi^. We find 
sizeable renormalization effects from scalar and pseudoscalar four-quark operators in the radiative 
decays and at 0(a s ) in hadronic b decays. Constraints on the Wilson coefficients of an extended 
operator basis are worked out. Further, the ratios Rh = B{B — ► Hfx + /B(B — ► He + e~), 
for = ifW, X s , and their correlations with the -B s — > n + /J,~ decay are investigated. We show 
that the Standard Model prediction for these ratios defined with the same cut on the dilepton mass 
for electron and muon modes, Rh = 1 + 0(m^/m^), has a much smaller theoretical uncer- 
tainty (< 1%) than the one for the individual branching fractions. The present experimental limit 
Rk ^ 1.2 puts constraints on scalar and pseudoscalar couplings, which are similar to the ones 
from current data on B(B S — ► /i+ We find that new physics corrections to Rk* and Rx s can 
reach 13% and 10%, respectively. 
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1 Introduction 



Flavor-changing neutral currents (FCNCs) are forbidden in the Standard Model (SM) at tree level 
and arise only at one loop. Hence, they are sensitive to quantum corrections from heavy degrees of 
freedom at and above the electroweak scale. The rare decays b — > 57, 5 — > sg and b — > sl + l~ , where 
i = e or fj,, are such promising probes. Measurements of these processes are rapidly improving by the 
present generation of B experiments and in the not too distant future by the Tevatron and the LHC. 
The analysis of b — > s transitions can be systematically performed in terms of an effective low-energy 
theory with the Hamiltonian (see, e.g., Ref. [1]) 

AC 

n cS = -^v^v^c^o^ + c^o'^)}. (i.i) 

The operators of* in Eq. (1.1) include dipole couplings with a photon and a gluon and dilepton 
operators with vector and axial- vector, as well as with scalar and pseudoscalar Lorentz structures. 
They are given as 1 

7 = ^m h {sa IH/ P R b)F' u ', O s = -m b {s a a^T a ap P R b p )G a ^ , 
9 S 9s 

2 2 

o 9 = e - (s^p L b) {hn) , O 10 = e -{si,P L b) , 

9 s 9 s 

2 2 
°s = 77T^ P « & ) W> °p = ^— 2 {sP R b){hd). (1.2) 

The operators (9i_ 6 can be seen in Ref. [5]. The primed operators in Eq. (1.1) can be obtained 
from their unprimed counterparts by replacing P L <-> P R . In the SM as well as in models with 
minimal flavor violation (MFV) where flavor violation is entirely ruled by the CKM matrix, the 
Wilson coefficients C[ are suppressed by the strange quark Yukawa coupling 

C[ ~ ^Q. (1.3) 

m b 

Furthermore, the SM contributions to scalar and pseudoscalar operators due to neutral Higgs-boson 
exchange are tiny even for taus since 

SM m e m b 

C s,p ~ — — ■ ( L4 ) 

Thus, in the context of the SM only the operators O 7 _i matter for semileptonic and radiative b — > s 
transitions. 



'Our definition of Os,p is different from that of Refs. [2-A] (i.e., without the factor of m^) in order for Cs,p to be 
dimensionless. As a consequence, the scalar and pseudoscalar operators have a non-vanishing anomalous dimension. 
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Our plan is to determine the coefficients c\ from a fit to the data and thereby testing the SM [6]. 
At present the number of measured independent observables is not sufficient, so one currently has to 
simplify the program and deal with a restricted set of operators. In this work we analyze the decays 
B -> X s j, B -> XJ + £~, B -> K^H+i', B s -> n + n~ with the following assumptions: 

(i) The effects of right-handed currents can be neglected, i.e., C[ ~ 0. 

(ii) The Wilson coefficients of scalar and pseudoscalar operators are proportional to the lepton mass 
C s ,p oc m e such that the coupling to electrons is negligible. This is automatically fulfilled if 
C s ,p are generated by neutral Higgs-boson exchange, but not in general within SUSY models 
with broken i?-parity. 2 

(iii) There are no CP-violating phases from physics beyond the SM. 

Therefore we take into account the Wilson coefficients C 7 _i and C s ,p- Model-independent anal- 
yses of the decays b — > S7 and b — > sl + l~ in the framework of the SM operator basis with 7 _ 1O have 
been previously performed in Refs. [8-10]. Distributions with an extended basis including O s ,p were 
analyzed for B -> X s i + t in Refs. [7, 1 1] and for B -> K*t+i~ decays in Refs. [12, 13] to illustrate 
possible new physics effects. In these works, however, no correlations between the just-mentioned 
decay modes and B s — > decays have been considered. In Ref. [3] the decays B s — > and 
B — > fi + fi~ have been studied model-independently. It has been shown that the Wilson coeffi- 
cients Cs,p can be of 0(1) while respecting data on the B s — > branching fraction, and thus are 
comparable in size to the vector and axial-vector couplings. For a combined study of B s — > ji + ji~ 
and B — > X s fi + fi~ decays in the minimal supersymmetric standard model (MSSM), see [14]. 

We perform here a combined analysis of the B s — > branching ratio and the observables 



where g^ ax = (m B — m K (,)) 2 for B — > K^*H + t~ and g^ ax ml for the inclusive decay modes. We 
also examine the low dilepton invariant mass region of the inclusive decays below the J/ip mass with 
<?max = 6 GeV 2 . Note that we use the lower cut of 4m 2 for both electron and muon modes in order to 
remove phase space effects in the ratio R H . Within the SM, we obtain clean predictions even for the 
exclusive decays 



2 Some i?-parity-violating SUSY models with horizontal flavor symmetries do have Cs.p oc mi. They can generate in 
general also helicity-flipped coefficients C- [7]. 




(1.5) 



RT = l + 0(ml/mi) 



(1.6) 
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Table 1: Branching fractions for various rare B decays [16-19]. The inclusive measurements as well 
as the corresponding theoretical predictions have been obtained for m e + e - > 0.2 GeV. The SM 
predictions are taken from Ref. [9] updated with B(B — > X c £u e ) = 0.108. 

Decay modes SM Belle BaBar 



X s e+e- (4.3 ± 0.7) x lO" 6 (5.0 ± 2.3tJ;?) x lO" 6 (6.6 ± 1.9±1;|) x lO" 6 

X s /i+/i~ (4-3 ± 0.7) x 10- 6 (7.9 ± 2.llfJ) x lO" 6 (5.7 ± 2&t\l) x lO" 6 

Ke + e- (3.6 ± 1.2) x 10~ 7 (4.8l^ ± 0.3 ± 0.1) x lO" 7 (7.4l^ ± 0.5) x lO" 7 

Kfj+fi- (3.6 ± 1.2) x 10~ 7 (4.81^ ± 0.3 ± 0.2) x 10- 7 (4.5™ ± 0.4) x lO" 7 

K*e + e~ (16.4 ± 5.1) x 10- 7 (14.91^1^ ± 0.2) x 10~ 7 (9.8^ ± 1.1) x 10~ 7 

B^K*fi + fi~ (12.4 ± 4.0) x 10" 7 (11.7±l;f ± 0.9 ± 0.5) x 10" 7 (12.7l£? ± 1.6) x 10 -7 



B 
B 
B 
B 
B 



which holds also outside the SM if C s ,p — 0. The normalization to the e + e _ mode in Eq. (1.5) was 
also discussed in Ref. [15] for the inclusive decays. 

This paper is organized as follows. In Sec. 2 we summarize the current experimental status and 
constraints on the decay modes of interest. Section 3 contains a discussion of new physics contri- 
butions to scalar and pseudoscalar four-quark operators and their impact on the Wilson coefficients 
of the SM operator basis. We investigate new physics effects in the decays b — > S7 and b — > sg. 
Model-independent constraints on the coefficients of the operators O7-10 in the presence of Os and 
Op are derived in Sec. 4. In Sec. 5 we study correlations between the branching ratios of the decays 
B s — > B — > X s £ + £~ and B — > K^*H + £~. In particular, quantitative predictions are obtained 

for the ratios R K (*) Xs - We summarize and conclude in Sec. 6. The anomalous dimensions, decay 
distributions for b — > s£ + £~ processes and auxiliary coefficients are given in Appendices A-D. 



2 Experimental status of b — ► s transitions 

We summarize recent results on the inclusive and exclusive b — > s£ + £~ decay modes in Table 1. 
These measurements are in agreement with the SM prediction [9] within errors. The experimental 
constraints we use in our numerical calculations are given below. Note that throughout this work we 
do not distinguish between B and B. 

(i) The combined results of Belle [19] and BaBar [16] for the inclusive b — > s£ + £~ decays yield 
the 90% confidence level intervals 

2.8 x 10~ 6 ^ B(B -> X s e + e-) < 8.8 x 10~ 6 , (2.1) 
3.5 x 10~ 6 ^ B(B -> X sfJl + fj,-) <: 10.4 x 10~ 6 . (2.2) 
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The statistical significance of the Belle (BaBar) measurements of B(B — > X s e + e ) and B(B — > 
X s fx + fi,~) is 3.4ct (4. Oa) and 4.7a (2.2a), respectively. To be conservative, we also use in our analysis 
the 90% C.L. limits [20] 

B(B -> X s e + e-) < 10.1 x 10~ 6 , (2.3) 
B(B -> X s n + n~) < 19.1 x 10~ 6 (2.4) 

and compare their implications with those of Eqs. (2.1) and (2.2). 

(ii) For the exclusive decay channels [17, 18] we obtain the following 90% C.L. ranges 

3.9 x 10~ 7 ^ B(B -> Ke + e-) ^ 7.7 x 10" 7 , (2.5) 

3.0 x 10" 7 ^ B(B -> Kfi+fT) < 6.5 x 10" 7 , (2.6) 
and 

6.5 x 10~ 7 ^ B(B -> K*e + e') < 17.9 x 10~ 7 , (2.7) 

6.7 x 10~ 7 ^ 5(5 -> K*n + n~) ^ 17.0 x 10~ 7 . (2.8) 

(iii) Using the experimental results displayed in Table 1 we find for the ratios Rh 

R Xs = 1.20 ± 0.52, = 0.81 ± 0.24, R K *\ no cut = 0.98 ± 0.38, (2.9) 

which translates into the 90% C.L. intervals 

0.34 < R Xs <: 2.06, 0.42 ^ 1.20, 0.35 ^ R K *\ QO cut ^ 1.60. (2.10) 

Here, i?x*| no cut is defined as with the lower integration boundary in the electron mode taken 
to be Ami, since experimental data on the B — > K*e + e~ branching ratios are published only for the 
full phase space region. We do not include effects from the small difference between the lower cut 
m e + e - = 0.2 GeV of the experimental analysis [16, 19] and 2m il used here. Furthermore, we neglect 
contributions to B — > Ke + e~ from the region below q 2 = Am 2 , where the rate is tiny due to the 
absence of the photon pole. The above ratios should be compared with the predictions of the SM 

flg* = 0.987 ±0.006, Rf 1 = 1 ±0.0001, R S ^\ no cut = 0.73 ± 0.01, (2.11) 

and 

i^| loW(?2 = 0.977 ±0.009, Rft = 0.991 ± 0.002, (2.12) 
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where "low g 2 " denotes a cut below 6 GeV 2 . The errors on the inclusive and exclusive ratios are due 
to a variation of the renormalization scale and of the form factors, respectively, see Sees. 4 and 5. 

(iv) The current world average of the inclusive b — > S7 branching ratio is [21, 22] 

B(B -> X sl ) = (3.34 ± 0.38) x 1(T 4 (2.13) 

with a photon energy cut E 1 > m b /20. 

(v) For the purely leptonic decays only upper limits exist. The branching ratio of the B s — > 
decay is constrained at 90% C.L. as [23] 

B(B S -> //+//") < 2.0 x 10~ 6 . (2.14) 

Note that there are preliminary 90% confidence level limits of 9.5 x 10~ 7 and 16 x 10~ 7 from CDF 
and D0, respectively [24]. 



3 New physics contributions to four-quark operators 

In this section we address the question whether new physics contributions to four-quark operators can 
spoil our model-independent analysis. Firstly, the QCD penguins appear in the SM and many 
extensions to lowest order only through operator mixing. They enter the matrix element of b — > 57 
and b — > si + i~ decays at the loop level. Hence, their impact is subdominant and new physics effects 
in QCD penguins are negligible for our analysis within current precision. Secondly, and this will be 
the important effect discussed in the remainder of this section, it is conceivable that the dynamics 
which generates large couplings to dileptons, i.e., to the operators Os,p, induces contributions to 
4-Fermi operators with diquarks as well. We introduce the following fermion / dependent operators 

O f L = (sP R b)(fP L f), O f R = (sP R b)(fP R f), (3.1) 

where for muons we identify the coefficients C£ R — e 2 / (16tc 2 )(Cs T Cp). We generalize here 
our assumption (ii) in the sense that the coupling strength is proportional to the fermion mass m/, 
which naturally arises in models with Higgs-boson exchange. In particular, the corresponding Wilson 
coefficients for b quarks proportional to m b can be potentially large. As will be discussed in the next 
section, current experimental data on the branching fraction of B s — > imply 3 

J\C b L (m w )\ 2 + \C R (m w ) | 2 ^ f mb(mw) y/2(\C s (m w )\ 2 + \C P (m w ) | 2 ) < 0.06. (3.2) 

Here, we anticipated our result in Eq. (4.4), i.e., an upper bound on \Cs,p\ and evolved according to 

d(C s ,p(fi) /m b (fj,))dfj, = 0, with the running 6-quark mass in the MS scheme given in Eq. (A. 2). 

3 In the MSSM with large tan j3 there are corrections to the down-type Yukawa coupling (see e.g. Refs. [25,26]). These 
corrections can be substantial in B decays, and have the form 1/(1 + et, tan (3) with |e&| < 0.01 [26]. 
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Figure 1 : Diagrams with an insertion of four-quark operators which contribute to the renormalization 
and the matrix element of the operator Og, and with an on-shell photon and no leptons to 7 . 

The Wilson coefficients C b LR are non-zero to lowest order interactions at the electroweak scale 
/j, ~ m w and can be significantly larger than the ones of the QCD penguins C 3 _ 6 (m 6 ) ~ 0(1CT 2 ). 
Hence, we have to study the potential impact of the operators O f LR on our analysis of b — > S7 and 
b — > sl + l~ decays. 

3.1 One-loop mixing with pseudoscalar and scalar operators 

Scalar and pseudoscalar four-quark operators enter radiative and semileptonic rare b — > s decays at 
one-loop level as shown in Fig. 1. To estimate their impact, we insert O b L and O b R into the penguin 
diagrams with an internal b quark and use fully anticommuting 75. The contributions from the diagram 
with closed fermion loop vanish by Dirac trace and by gauge invariance or vector current conservation, 
i.e., after contraction with the lepton current. For simplicity, we work in the "standard" operator basis 
Oi given in Appendix A. We obtain non- vanishing contributions from O b R and O h L to 7 and Og, 
respectively. The diagrams with an internal s quark contribute to the helicity-flipped coefficients. 
They are suppressed by a factor m s /m b and therefore can be neglected. We obtain the following 
corrections to the Wilson coefficients at the scale /i b = m b 

5C 7 (m b ) = ];\n^fC b R (m w ), (3.3) 
6 ml 

5Cg(m b ) = l\n^C b L (m w ). (3.4) 
9 m h 

These infinite renormalization contributions survive in the limit a s — > 0, which is similar to what 
happens in the SM for the mixing of 2 onto Og [27]. With the upper bound in Eq. (3.2) we find that 
the new physics effect from O b LR is small, of the order of one percent for Og, but a few x 0(10%) for 
O7. The reason is simply that Cj M (m b ) is more than an order of magnitude smaller than Cg M (m b ), 
which in addition has a smaller anomalous dimension. 
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Other operators contributing in the SM but subleading in the decays b — > 57 and b — > s£ + £~ are 
also subject to similar new physics effects. To be specific, the Wilson coefficients of the chromomag- 
netic dipole operator and the QCD penguin operators receive corrections from the diagrams in Fig. 1 
with diquarks instead of leptons and the intermediate photon replaced by a gluon. We find 

SC 8 (m b ) = -^m^f C b R (m w ), (3.5) 

SC 3 , 5 (m b ) = -1^ In ^C b L (m w ), (3.6) 

6C,, 6 (m b ) = In ^Ci(m^), (3.7) 
6 An m b 

which are relevant to hadronic B decays. 4 Quantitatively, the renormalization of the gluon dipole 
operator can be order one. (We study this in more detail below.) The impact on the QCD penguins 
can be up to several percent. As mentioned earlier, new physics contributions to the operators 
are subdominant in b — > s-y and b — > s£ + £~ decays. Since the renormalization of Og by scalar and 
pseudoscalar operators is small, too, we can safely neglect the effects of induced four-quark operators 
of the type O h L in our analysis of semileptonic and radiative b — > s decays. We remark that scalar and 
pseudoscalar operators also mix with the electroweak penguin operators Oj_ lQ (see Appendix A) at 
order a/ An. We have calculated for completeness the corresponding anomalous dimensions, which 
can be seen in Appendix B. 

To get a more accurate estimate of the new physics corrections to the magnetic penguin opera- 
tors, we resum the leading logarithms in Eqs. (3.3) and (3.5) by means of the renormalization group 
equations in the MS scheme [1]. Both operators O b R and O h L induce additional operators under renor- 
malization (see Appendix B). The anomalous dimensions of each set are known at next-to-leading 
order (NLO) [29], with no mixing between the sets. We have calculated the leading-order mixing of 
O h L R onto 03-9. 5 The anomalous dimensions are given in Appendix B together with the respective 
leading-order self-mixing of both O b R and O h L sets. Numerically, we obtain 

5C 7 (m b ) ~ 0.7lC b R (m w ), 5C s (m b ) ~ -2.95C b R (m w ), (3.8) 

which implies sizeable contributions to the branching ratios of the radiative decays. We study the 
phenomenology in Sec. 3.2. 

4 The decay B — > <f>Ks has been studied in Ref. [28] including 0(a s ) corrections to the matrix element. The leading 
logarithmic contributions in Eqs. (3.5)-(3.7), however, have not been taking into account, which explains the huge fi 
dependence found in these papers. We checked that the ln(m{,//i) terms of the 0(a s ) corrections are canceled by the 
contributions in Eqs. (3.5)-(3.7). 

5 The computation of the anomalous dimensions at NLO is being performed in Ref. [30]. 
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The mixing of scalar and pseudoscalar operators in Eq. (B.l) onto the dipole operators has been 
studied previously in the context of the two-Higgs-doublet model [31] and in supersymmetry with 
gluino contributions to b — > S7 [32]. While our results agree with the ones presented in Ref. [32], 
they are at variance with those given in Ref. [31]. In particular, we disagree with the conclusion made 
therein that the scalar and pseudoscalar operators do not mix with Og. 

3.2 Implications for the decays b — > sj and b — > sg 

We now investigate the phenomenological consequences of the mixing effects presented above for 
radiative B decays. To illustrate how large these corrections can be, we normalize the Wilson coef- 
ficients in the presence of new physics to the ones in the SM, and denote this ratio by £, such that 
£ SM = 1. We obtain to next- to-leading order in the SM operator basis and to leading logarithmic 
approximation in C b R 

£ 7 {m b ) = 0.514 + 0.450 fr(m w ) + 0.035 &(m w ) - 2.319 C b R (m w ), (3.9) 

£ 8 (m b ) = 0.542 + 0.458£ 8 (m w ) + 19.790 C b R (m w ). (3.10) 

Given the upper bound in Eq. (3.2) corrections of up to 14% and 119% to £ 7 and £ 8 can arise. We work 
out correlations between £ 7 and £ 8 from B(B — > X s ~f) given in Eq. (2.13) and B(B — > X s g) < 9% 
at 90% C.L. [33], using the analytical formulae of Refs. [34-36]. We obtain the allowed regions at 
the /j, b scale shown in Fig. 2 for C R (m w ) = (left plot) and C b R (m w ) = 0.06 (right plot). The 
theoretical uncertainty from the prescription of the charm-quark mass has been taken into account by 
including both solutions obtained for m c /m b = 0.22 and 0.29 [37]. From Fig. 2 we see that A 7 = 
for C R (m w ) = 0.06 is allowed by present data on the b — > sg branching fraction. This particular 
scenario could be excluded by an improved experimental analysis of b — > sg. Also, if C R (mw) is 
near its upper bound, it implies a contribution to the matching conditions for Ci^{mw) in order to be 
consistent with experimental data. 

In summary, we find that the impact of C\ on the matrix element of Og is small, at most a few 
percent, and thus can be neglected. On the other hand, contributions to the dipole operators are in 
general non-negligible. They can be avoided assuming C R (m w ) ~ 0, i.e., 

C s + C P = 0. (3.11) 

In the remainder of this work we discuss the phenomenology with and without this constraint. Note 
that the absence of logarithms in the matching conditions for C 7iS (mw) from neutral Higgs-boson 
exchange in a two-Higgs-doublet model type II [26, 38] is consistent with the fact that in this model 
Eq. (3.11) is satisfied [2]. This is also the case for the MSSM with MFV at large tan (3 [3]. 
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Figure 2: Constraints on £ 7fi (jjt b ) from £(& -> s 7 ) for C£(M W ) = (left plot) and C£_(M W ) = 0.06 
(right plot). We also show the upper and lower bounds on £s (/•*&) f° r me experimental limit B(B — > 
X s g) < 9% [33] (dashed lines) and for an assumed value of B(B — > X s g) < 3% (dash-dotted lines). 

4 Model-independent analysis 

In this section we give the theoretical framework that we use to analyze the decays B — > X s ^, 

— > . We then work out model-independent constraints on 
the coefficients of the operators C7-10 and C^p. 

4.1 Wilson coefficients and matrix elements 

The matrix element of inclusive b — > s£ + £~ decays contains contributions from the photon dipole 
operator O7, the dilepton operators Cg^o and in models beyond the SM also from Os,p- The decay 
distributions in the SM are known to next-to-next-to-leading order (NNLO) [5,39-41], which corre- 
sponds to NLO in b — > 57. We use the NNLO expressions for the operators 07,9,10 and lowest order 
ones for Os,p since 0(a s ) corrections to the matrix elements of leptonic scalar and pseudoscalar 
operators in these decays are not known. 

Further, we assume that the contribution from intermediate charmonia has been removed with 
experimental cuts. Non-perturbative corrections [42] affect the branching ratio by at most few percent 
and we do not consider them here. We neglect the mass of the strange quark but keep the muon mass 
consistently, because according to our assumption (ii) also C s ,p counts as one power of mi and can 
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be enhanced in models beyond the SM. 

The dilepton invariant mass spectra for inclusive and exclusive b — > s£ + £~ decays are given in 
Appendix C. The effective coefficients which enter the decay distributions are written as [5, 39] 



C 7 e 



cff 



1 H ^7{S) 



TT 



AM ~ 



47T 



i=l,2 



(7)f 



S)C7<%) + ^''(S)^) 



(4.1) 



CefF 
9 - 



1 H a; 9 (s) 



7T 



[AM + T 9 h(m 2 c , s) + U 9 h(l, s) + W 9 h(0, s)] 



An 



=1,2 



°10 



1 H w 9 (s) 



7T 



(4.2) 



(4.3) 



where m c = m c /m b , s = q 2 jm\ and A i: T 9 , C/ 9 , Wg are given in Appendix D. The function h(z, s) 
originates from the one-loop matrix elements of the four-quark operators 0i_ 6 (see Fig. 1) and can be 
found in Ref. [5]. The functions Lu i: arise from real and virtual a s corrections. They can be seen 
in Refs. [5,39] together with uj 79 which replaces uo 7 and lu$ in the interference term Re (C^Cf 1 *) in 
the decay rate. In the calculation of the decay rate we expand in powers of a s and retain only linear 
terms. Note that the uoi include only that part from real gluon emission which is required to cancel the 
divergence from the virtual corrections to the matrix element of the C^. Further gluon bremsstrahlung 
corrections in b — > sl + l~ decays [40,43] are subdominant over the whole phase space except for very 
low dilepton mass and are not taken into account here. In our numerical analysis we choose a low 
value for the renormalization scale, fj,b = 2.5 GeV, because this approximates the full NNLO dilepton 
spectrum by the partial one, i.e., with the virtual 0(a s ) corrections = in Eqs. (4.1) and (4.2) 

switched off [9]. This is beneficial since the are known in a compact analytical form only for the 
low dilepton invariant mass region [39]. For the exclusive B — > K^i + i~ decays we set uo,i = 0, 
since these corrections are already included in the corresponding form factors. We do not take into 
account hard spectator interactions [44]. 

Below we work out model-independent bounds on = Af M + Af p . They differ from the "true" 
Wilson coefficients C, by penguin contributions that restore the renormalization scheme independence 
of the matrix element [27]. In addition A 9 contains logarithms from insertions of the four-quark 
operators Oi- 6 into the diagrams of Fig. 1. Explicit formulae relating A^ and are given in Appendix 
D. As discussed in Sec. 3, we neglect new physics contributions to the QCD penguin operators. 
In our numerical study we use fs a = 200 MeV and 238 MeV [45] and the parameters given in 
Table II of Ref. [9] except for B(B -> XJu e ) = 10.80% [46]. Form factors and their variation 
are taken from Ref. [10]. We give the SM values for completeness: A^ M (2.5 GeV) = -0.330, 
Al M {2.5 GeV) = 4.069 and Af™ = -4.213. 
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4.2 Constraints from B s —> ii + fi 

An upper limit on the B s — > /i + 'fi~ branching ratio constrains the scalar and pseudoscalar couplings 

VICsW + ICpM + M/OI 2 < 3.3 \ bUJ ' r 11 '' 



x 



\V tb Vt*s 
0.04 



2.0 x 10- 6 

2 



4.4 GeV 



238 MeV 



1/133 



(4.4) 



Here, we neglected the factor (1 — Am^/rn^J in front of |Cs| 2 , see Eq. (C.3), and defined 5i (fi) = 
2m fJi m b (fx) /m 2 Bs A 10 . The bound given in Eq. (2.14) also implies the upper limits 

B(B S -> e + e~) ^ 4.7 x 10 -11 , B(B S -> r+r") < 4.2 x 10" 4 . (4.5) 



4.3 Constraints from b 



S'J 



The measured — > 57 branching fraction puts constraints on the dipole operators. In the absence 
of scalar and pseudoscalar couplings C b R (see Sec. 3), which renormalize both electromagnetic and 
gluonic operators, the two solutions A 7 (fj, b ) ~ ±Aj M (fi b ) are allowed. This is the case if Eq. (3.11) is 
satisfied. We update the NLO analyses of [9, 35] with the inclusive b — > S7 measurement in Eq. (2. 13) 
and B(B -> X c eu e ) = 10.80% and obtain the ranges (fi b = 2.5 GeV) 

-0.36 < A 7 < -0.17 or 0.21 ^ A 7 ^ 0.42. (4.6) 

The corresponding correlation between A 7 and A 8 can be seen in the left plot of Fig. 2. For C b R (m w ) = 
0.06, on the other hand, the experimental constraints on A 7 are much weaker (right plot of Fig. 2). 



4.4 Constraints from b 



In the presence of new physics contributions proportional to the lepton mass we use data on the 
electron modes to constrain the dilepton couplings A 9;10 . From the upper bound on B(B — > X s e + e~) 
given in Eq. (2.3) we obtain 

' 9.0 for A 7 < 
8.9 for A 7 > 



A 



-0.58 1 2 
9+1.05| 



+ \A l0 \ 2 ^ 



y / |A 9 + 0.15| 2 + |A 10 | 2 ^ 9.1 for A 7 = 0. 
The range on B(B — > X s e + e~) given in Eq. (2.1) yields upper and lower bounds 

U for A 7 < 



(4.7) 




a —0.58 I 2 
^9+1.05 I 



Aol 2 ^ 



8.3 for A 7 > 
4.8 < y / |A 9 + 0.15| 2 + |A 10 | 2 < 8.5 for A 7 = 0. 



(4.8) 
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Similar bounds can be obtained from data on the muon modes together with the upper limit on Cs,p 
in Eq. (4.4). The lower limit on B[B — > X s [i + [i~) in Eq. (2.2) yields 

4 -i.4|2 , i 4 / 3.8 (3.5) for A 7 < 

A9+1 - 9 ' ^ \ 3.5 (3.2) for A 7 >0 

|A 9 + 0.15| 2 + |^io| 2 > 4.7(4.4) for A 7 = (4.9) 

for f Bs = 238 MeV (200 MeV). Our constraints on A 9A0 given in Eqs. (4.7)-(4.9) are displayed in 
Fig. 3. Like in the analysis with the restricted SM basis in [9], A 9 = A w = is excluded even in the 
presence of new scalar and pseudoscalar interactions. 



4.5 Constraints from R K 

The experimental bound R K ^ 1.2 in Eq. (2.10) provides constraints on the scalar and pseudoscalar 
Wilson coefficients complementary to those from the B s — > branching fraction given in 

Eq. (4.4). Varying A 7i9jl0 according to Eqs. (4.6), (4.7)-(4.9) we obtain (fj, b = 2.5 GeV) 

\/|Cs| 2 +|Cp + A 10 | 2 <c4.5. (4.10) 

Here, Ai stems from the interference term of Cp and A w in the B — > Kfi + fi" rate, see Eq. (C.4), 
which can be neglected for large values of Cs,p- If the bound on R K improves e.g. to 1.1, then the 
value on the r.h.s. of the above equation changes to 3.2. 



5 Correlation between B s — ► and b — ► s£ + £ decays 

In this section we study correlations between the ratios Rh defined in Eq. (1.5) and B(B S — > ji + 
We restrict ourselves to the case Cs = —Cp, hence a vanishing A 7 is excluded as shown in Sec. 3.2. 
We further assume that A 9)10 are SM valued while A 7 is allowed to vary in the intervals given in 
Eq. (4.6). This particular scenario is, for example, realized in the MSSM with MFV at large tan (3. 
The maximum values of Rh are summarized in Table 2 of Sec. 6 for different new physics scenarios. 

The correlations depend sensitively on the decay constant of the B s meson. We display our results 
fo r fB s — 200 MeV and 238 MeV except for the inclusive decays, where we vary between these two 
values. As described in Sec. 4 we use the partial NNLO expressions. Therefore, the plots are obtained 
for fixed renormalization scale fi b = 2.5 GeV. For the analysis of the exclusive decays we show the 
uncertainty from the form factors. 

The SM predictions for the ratios R H are summarized in Eqs. (2.11) and (2.12). The theoreti- 
cal uncertainty for the inclusive decays is due to the variation of the renormalization scale between 
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-10 



-10 



A 7 = 




-5 





A 9 (2.5 GeV) 



10 



Figure 3: Allowed regions in the Aq-A w plane in the presence of scalar and pseudoscalar operators 
from data on inclusive b — ► s£ + £~ and b — > S7 decays for different values of A 7 . The shaded areas 
are obtained from the upper bound on B(B — > X s e + e~) and the lower bound on B(B — > X s n + fi~), 
Eqs. (4.7) and (4.9) with fs a = 200 MeV. The two remaining contours indicate the allowed regions 
from the 90% C.L. measurement of B[B — > X s e + e~) given in Eq. (4.8). Since the bounds from 
B(B — > X s fi + fi~) for f Bs = 238 MeV give very similar results, we do not show the corresponding 
contours. 
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2.5 GeV and 10 GeV. Since we are using the partial NNLO expressions this small error below one 
percent on R^ M might even be overestimated. For comparison, we give the corresponding numbers 
at NLO R s ™' NLO = 0.974 ± 0.006 and i^' NL0 |i w ? = 0.972 ± 0.005. The SM prediction for the 
B s — > branching ratio is (3.6 ± 1.4) x 10 -9 , where the main theoretical uncertainty results from 

the B s decay constant. It can be considerably reduced once the B®-B® mass difference is known [47]. 



5.1 Exclusive B — > K£ + £ decays 



Figure 4 shows the correlation between R K and the B s — > fi + ji~ branching ratio for two values of 
the B s -meson decay constant and different signs of A 7 and Cp. As illustrated by the solid lines in the 
upper left plot, the dependence of Rk on the form factors is very small and hence this observable is 
useful for testing the SM. For comparison, the uncertainty on the B — > K£ + £~ branching fraction due 
to the form factors is ~ 30% [9]. While being consistent with B s — > data given in Eq. (2.14), 

an enhancement of R^ 1 by ~ 60% is excluded by the current upper limit on Rk (dotted lines). 

Furthermore, the ratio Rk provides a bound on C$,p which is competitive with the limit from 
B(B S — > yu + /i~) in Eq. (4.4). For two values of Rk we find {ji b = 2.5 GeV) 



^ 2.3 for R K = 1.1, 



whereas data on B s — > decays give 

V\Cs\ 2 + |C P -0.15| 2 < 3. 



(5.1) 





1/2 


"238 MeV" 


2.0 x 10" 6 




fs 3 



(5.2) 



We recall that R K = 1.2 corresponds to the current 90% C.L. upper limit, see Eq. (2.10). 



5.2 Exclusive B -> K*i + i~ decays 

The results for R K * versus the branching ratio of B s — > are shown in Fig. 5. Note that the 

variation from the form factors is much larger than in Rk- This is caused by the form factor A , which 
drives the Cs,p contributions to Rk*- Its theoretical uncertainty in light cone QCD sum rules [10], 
which we use in our analysis, is twice as large as in f relevant for R K . New physics effects in R K * 
can be as large as 30% [allowed by B s — > data in Eq. (2.14)] but are restricted to be less than 

12% once data on R K are taken into account. For the ratio with no lower cut on the electron mode we 
find including all constraints R K * | no cut ^ 1.01, an enhancement of 36% over its SM value. 
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A 7 < 






C P < 










^^^T 238 


RSM 





0.5 



1.5 

10 6 




1 1.5 

->• //+//-) x 10 6 





1 1.5 2 

B(B S -> x 10 6 

Figure 4: Correlation between i?^ and the B s — ► branching ratio for different signs of and 

Cp, two values of fs s in MeV and Ag io = ^4gYo- The shaded areas have been obtained by varying 
the B — > K form factors according to Ref. [10] and A 7 as given in Eq. (4.6). In the upper left plot, 
the form factor uncertainty is illustrated for fixed A 7 = A 7 M and Jb s = 200 MeV by solid lines. The 
dotted lines correspond to the 90% C.L. upper limit on Rk in Eq. (2.10). Dashed lines denote the SM 
prediction for R K . 
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09 0.5 " 1 " 1.5 2 09 0.5 1 1.5 2 

B(B a -S- /J+/J,-) x 10 6 B{B a -S- x 10 6 

Figure 5: Correlation between i?^* and the B s — > branching ratio (see Fig. 4 for details). 

The dotted lines represent the maximum value of Rk* consistent with the experimental upper limit 
Rk ^ 1-2. Dashed lines denote the SM prediction for Rk*- 
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5.3 Inclusive B — > X s l + l decays 

In Fig. 6 we show the correlation of R Xs with the B s — > branching ratio for the full spectrum 

with s max ~ 1 (upper plots) and for the low dilepton mass with s max = 0.26 (lower plots). Order one 
effects in Rx s from scalar and pseudoscalar interactions are excluded by current data on B s — > 
contrary to the results of Ref. [15] but in agreement with Ref. [48]. We find a maximum value of Rx 3 
of 1.08 (full spectrum) and 1.05 (low dilepton mass) from the experimental upper limit on Rk- These 
bounds on the B — > X s /j, + fx~ branching ratio are similar to the ones from B s — > data previously 

obtained in [48]. While an enhancement of the B — > X s /j, + fx~ branching ratio of 0(10%) is within the 
uncertainty of the SM prediction, a corresponding effect in the ratios R Xs can be well distinguished 
from the SM ones. 

6 Summary 

We performed for the first time a model-independent analysis of b — > s processes in an extended 
operator basis, the SM one with O 7 ^i plus scalar and pseudoscalar operators Os,p with dileptons. In 
our phenomenological analysis we took into account experimental constraints from inclusive b — > 57, 
b — > s£ + £~ (£ = e, //) and B s — > fi + ji~ decays. Further, we used data on the ratio of B — > Kfi + fi~ to 
B — > Ke + e~ branching ratios, Rk- We made a few assumptions to facilitate this analysis: no right- 
handed currents, the couplings to the scalar and pseudoscalar operators are driven by the respective 
fermion mass and no CP violation beyond the CKM matrix. 

We studied the effects of scalar and pseudoscalar operators involving b quarks O h L R . Already at 
zeroth order in the strong coupling constant these operators mix onto the SM basis: O h L proportional 
to Cs — Cp onto the 4-Fermi operators with dileptons and O b R proportional to Cs + Cp onto the 
photonic and gluonic dipole operators. Furthermore, we find that the QCD penguins get renormalized 
at 0(a s ) by O h L . While being negligible in b — > s£ + £~, these corrections are important for hadronic 
b decays. In particular, they cancel the strong /j, dependence of the B — > <j>K s amplitude reported 
recently in Ref. [28]. The lowest order anomalous dimensions involving O b R and the dipole operators 
have been calculated before in Ref. [32], whereas the ones with O b L and the 4-Fermi operators are a 
new result of this work. Numerically, we find that for Cs = —Cp the effects of O b L R are negligible 
for our model-independent analysis. However, for C s 7^ —C P there is a significant impact from scalar 
and pseudoscalar couplings on the dipole operators. In particular, the branching ratio for the decay 
b — > S7 can be obtained completely without any contribution from the electromagnetic dipole operator 
O7. This rather extreme scenario could be excluded by improved data on the b — > sg branching ratio, 
as illustrated in Fig. 2. Except for the case of A 7 ~ 0, the bounds we obtain on the coefficients A 910 
are similar to previous results in the SM operator basis [9]. The non-trivial renormalization effects we 
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0.5 1 1.5 2 0.5 1 1.5 2 

B{B, -s- n+n-) x 10 6 B(B S -> x 10° 

Figure 6: The dependence of Rx a on the 5 S — > branching ratio for different signs of A 7 and 

C P , and Ag io — Aff . The upper plots correspond to the whole dilepton invariant mass spectrum 
while the lower ones correspond to the low-g 2 region as described in the text. The shaded areas have 
been obtained by varying f Bs between 200 MeV and 238 MeV and A 7 according to Eq. (4.6). In the 
left plots, the solid lines indicate the uncertainty from the variation of fs s for fixed A 7 = Aj M and 
Cp < 0. The dotted lines represent the maximum allowed values of Rx s obtained for Rk ^ 1.2. 
Dashed lines denote the SM prediction for Rx a - 
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Table 2: Upper bounds on the ratios Rh for different Cs,p scenarios for A 9 ^ being SM-like and in 
parentheses without this constraint. Data on b — > S7, 6 — * s£ + £~, i? s — > and i?^ are taken into 

account. 



Ratio 


SM 


C s ,p = 


Cs — —Cp 


C s ^ —Cp 


Rk 


1.00 


1.00(1.00) 


1.20 (1.20) 


1.20 (1.20) 


Rk* 


0.99 


1.00(1.00) 


1.11 (1.12) 


1.12(1.12) 


R K*\no cut 


0.74 


0.91 (0.97) 


1.01 (1.07) 


1.11 (1.12) 


Rx s 


0.98 


0.99 (0.99) 


1.08 (1.08) 


1.08 (1.08) 


low g 2 


0.97 


0.99 (0.99) 


1.05(1.06) 


1.07(1.07) 



encountered show that a model-independent analysis can be quite involved in an enlarged operator 
basis. 

We worked out correlations between the ratios Rh defined in Eq. (1.5) and the B s — > ji + pT 
branching ratio for Cs = —Cp and A 9 W being SM-like, summarized in Figs. 4-6. This particular 
scenario also applies to the MSSM with MFV at large tan/3. Figure 4 shows that a bound on Rk 
implies a bound on B(B S — > /x + /U~) and vice versa. Current data on these observables yield very 
similar constraints on Cs,p given in sections 4 and 5. Note that in the above-mentioned MSSM 
scenario B(B S — > and B° s -B° s mixing are correlated [26]. A similar correlation between R K 

and in general with larger theoretical errors also with the other i?'s and B®-B® mixing then holds 
in this model, too. We stress that in our analysis we take into account information on branching 
ratios only from inclusive decays. The data on exclusive decays enter our analysis only via R K 
which depends only weakly on the form factors, as can be seen from Fig. 4. The largest theoretical 
uncertainty in the correlations is due to the £? s -meson decay constant. 

We further calculated the maximal allowed values of the ratios R H , summarized in Table 2. Since 
we use the partial NNLO expressions for the Wilson coefficients, they have been obtained at the 
scale ii b = 2.5 GeV. We see that large, order one corrections to the respective SM values are already 
excluded. Note that these upper bounds are insensitive to fp a because current data on R K ^1.2 are 
here more constraining than B(B S — > u, + fi~). The effect from C s ,p on B — > K decays is always 
bigger than on B — > K* and B — > X s decays. The reason is that besides different hadronic matrix 
elements in these decays the photon pole |A 7 | 2 /s, which is absent in the B — > K decay, dominates 
the rate for very low dilepton mass. The inclusive decay with the spectrum integrated only over the 
low dilepton invariant mass is even less sensitive, since the lepton-mass-dependent contributions are 
suppressed by small s, see Eq. (C.18). 

Contributions from scalar and pseudoscalar operators with V + A handedness can be included in 
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the B s — > branching ratio by Cs,p — > C^p — C5 P and into the Z? — > K£ + £~ spectrum by 
Cs,p Cs,p + Cgp- Hence, the correlations we presented between B(B S — > n + and i?x break 
down if both chirality contributions Cs,p and C^p are non- vanishing. Since Rk constrains the sum 
and B s — > the difference of the coefficients, combining these two [Eqs. (4.4) and (4.10)] yields 

an upper bound on the magnitude of the individual coefficients of \C$ P \ ^ 4.3. This excludes large 
cancellations and holds even with right-handed contributions to the SM operator basis. 

In conclusion, b — > si + i~ induced decays can have a splitting in the branching ratios depending 
on the final lepton flavor from physics beyond the SM. Hence, averaging of electron and muon data 
has to be done carefully in order not to yield a model-dependent result. The effect from scalar and 
pseudoscalar couplings is best isolated in the theoretically clean observables R H with the same cuts 
on the dilepton mass. On the other hand, the ratio R K * | no cut constructed with physical phase space 
boundaries is also sensitive to new physics not residing in Cs,p, as can be seen from Table 2. 

Note added. The lowest order mixing of scalar and pseudoscalar operators onto the SM basis 
calculated in Sec. 3.1 has been taken into account in a revised version of the first paper of Ref. [28]. 
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A Standard operator basis 

In this appendix we give the "standard" operator basis [1] 

Ox = (s a ^P L c p )(c^P L b a ), d 2 = (s all ,P L c a )(c^P L bp), 



3 = (s a ^P L b a ) (Qf>-f p LQfi), A = (s al »P L bfi) (Qp-f p LQ a ), 

q=u,d,s,c,b q=u,d,s,c,b 



o 5 = (s al ^p L b a ) &-f p Rqp)> o & = {s alll p L b p ) (Qtn li PRQ<*), 

q=u,d,s,c,b q=u,d,s,c,b 
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o 7 = 



m b (s a a^P R b a )F^, 8 = J^m b (s a a^T^P R bp)G a ^ , 



167T 2 

2 



167T 2 

-,2 



(s al ,P L b a )(!^£), O w = -^(s a ^P L b a )(£ri5t). 



(A.1) 



lOTT^ 1671^ 

Here Q g denotes the charge of the q quark in units of e, a, (3 are color indices, a labels the SU(3) 
generators, P LR = (1 =F 7s)/2 and m 6 = m b (fi) is the running mass in the MS scheme, 



m h (jj) = m 



with 7 <? = 8, >y£> = 1012/9, (3 = 23/3, (3 X = 116/3. 



pole 
b 



X 1 ) 



_ a s (m^ ole ) 16 
3 



(0) 
-, ™ 

2/3 



1 + 



7# At 



(0) 



2/9o 2/5 2 



Oi s (n) - a s {mf le ) 



47T 



(A.2) 



B New operators and mixing 

The new physics operators containing scalar, pseudoscalar and tensor interactions are written as 

<5ll = (s a P R b a )(b a P L b a ), 012 = {s a P R bp){bpP L b a ), 

0i3 = (s a P R b a )(b a P R b a ), 0i4 = {s a P R bp){b p P R b a ), 



015 = (s a a^P R b a )(b a( r^P R b a ), 0i 6 = (s^PMibpa^PM, (B.l) 

where a^ u = (i/2)[7 M , 7„] and 011,13 = 0^ ^ in Eq. (3.1). For completeness, we give their lowest 
order self mixing [29,32,49], i.e., among On, i2 



7 



a. 



■16 



4tt \ -6 2 

and among 13 ,...,0 16 

/-16 1/3 -1 \ 
a s -6 2 -1/2 -7/6 
4tt 16 -48 16/3 
V-24 -56 6 -38/3/ 



(B.2) 



7 = 



(B.3) 
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We obtain the following lowest order anomalous dimensions for the mixing of O13, . . . , 0±q onto 7:& 
713-iej = Qd(l, N c , -[4 + 8N C ], -[AN C + 8]), 7l3 _i 6)8 = (1, 0, -4, -8), (B.4) 
and ofOn.12 onto 3 ,...,d 6 ,0 9 

711,3-6 = -7^2 -1, ^r, -1 I , 712,3-6 = 0, 7n, 9 = — , 7i 2 ,9 = , (B.5) 

where iV c is the number of colors. Note that Eq. (B.4) is in agreement with [32]. For the mixing of 
On, 12 onto the electroweak penguins Oj, . . . , O\ , we find 



7xVio = -^^(1,0,1,0), 7l Vio = -£^^(1,0,1,0). 
The remaining leading order anomalous dimensions vanish. 



(B.6) 



C Differential decay distributions 

We neglect the s-quark mass and introduce the notation 



/ / AttV 

rhi = rrii/mB, s = q 2 /m B , u(s) = J A ( 1 

A = 1 + m A K{t) + s 2 - 2s - 2m 2 KM (l + s) 
for the exclusive decays and 

m t = m t /mf c , s = q 2 /(m p b ° le ) 2 



(C.l) 



(C.2) 



for the inclusive modes. Then, the various decay distributions in the presence of scalar and pseu- 
doscalar operators can be written as follows. 



C.l B, 



£ + t 



£2 2 3 f2 

r( 5s _> i+e-) = — - B - JBs \v tb v; s \ 2 



X 



64tt 3 
Am 2 

1 7T- 



m. 



m. 



m Bs C s 



m b 



m Ba C P 2m e 

1 A w 



(C.3) 



with A 10 defined in Eq. (D.4) and C s ,p = C s ,p{^), m b = m b (u). 
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C.2 B -► 



d r( S ^i f «-)_G| a W B|w , |2 . ( . ) . () ^, j2 + |c , |2 , 



ds 



2 10 7T 5 



A - 



MS 



+ A\C'\ 2 m 2 E (2 + 2m^ - s) + 8Re (C"D'*)m?(l - m\) + 4| J D'| 2 m^s + |T P | 2 s + \T s \ 2 (s - Am 2 e ) 
+ 4Re (D'T^)m e s + 4Re (C"7£)m/(1 - m 2 K ) 1, (C.4) 



with 



A' = Cf(S)f + (s) + 



2m 



1 + 



<T/r(S), 



C = C*U{§), 



^2 



, 1 — m 



s,p 



m b 



Cs, P fo(s), 



(C.5) 
(C.6) 
(C7) 

(C.8) 



where the definition of the form factors can be found in Ref. [10]. The Wilson coefficients Cf s can 
be obtained from the ones in Eqs. (4.1)-(4.3) with cj 7i g i7 g = 0. 

C.3 B -> K*£+£- 



dT{B -+ KHH-) _ G 2 F a 2 m\ 2 f 1 



\ A ^sx(l + ^) + \E\ 2 su(s) 2 



+ 



+ 



Am 2 K , 



A 



Am 2 K ~ M 
1 



\B\ 2 
\C? 



A - + 8m 2 K , (s + 2m|) 



+ \F\ 



A 



V s 
u(s) 2 



+ 8m 2 K *(s — Am 2 ) 



A- 



A - + 4m 2 (2 + 2m£, - s) 



2m^* 
+ 4m;; A 



Re(BC*) 



A- 



(1 -rh 2 K , - s) + Re (FG*) 



A 



A- 



(1 - m\, - s) 



A 



- 2m|[Re (FH*) - Re (GH*)(1 - rh 2 K ,)}^- + m 2 e \H\ 2 s^~ + \X P \ 2 s 



A 



m 



K* 



Am 2 K , 



+ |X 5 | 2 (s-4m 2 )-^-m^[Re(FX|,)-(l-m^0Re(G'X;)-sRe(M*)]-4-^ (C.9) 



Here, 

A = 



1 + m K 



' Am 2 K * 



-Cf(s)V(s) + ^Cf Tl (s), 



m" K * 



(CIO) 
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B = (l+m K *) 

c 1 



Cf(s)A,(i) + 5^(1 - m K .)Cf T 2 (s) 
s 



1 — m\, 



(1 - m K *)Cf{s)A 2 {s) + 2m 6 C? 



eff 



Us) + 1 -4^-T 2 {s) 



E 



1 + m,K* 
F=(l+m K *)C$A 1 (s), 



G 



1 + rh K * 



1 



H = ~ c To K 1 + %.)Ai(s) - (1 - m x *)^2(s) - 2m K *A (s)} , 
s 

X s ,p = - 2 -^A (s)C s , P , 
with the form factors defined in Ref. [10]. The (7f ff 's are given in Eqs. (4.1)-(4.3) with uo 7 $ 



C.4 B XJ+t 



dT(B -> X s £ + r) G%a 2 (m : 



2 P°le\5 



X 



ds 3 x 2 8 tt 5 

12Re (CfCf*) + 



IW(i-iNi 



4m| 



+ (|c 9 cf T + |CTol 2 ) 



1 + 2s + 



2mf (1 - s) 



3. 
+ 2 S 



1 - 



4m^ 



\Cs\ 2 +\C P \ 



+ 6m,Re (C P C^*) |>, 



with Cf defined in Eqs. (4.1)-(4.3). [Equation (C.18) agrees with Ref. [7] for m s = 0.] 



D Auxiliary coefficients 



a s {fi) 3 9 3 9 
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- - £^ CM + £ QOifri? In ^ + |c7 3 (/i) + y C 5 (/x) + ^C 6 (/z), (D.3) 
Air 

A 10 = —^CM, (DA) 

T 9 = ^d(Ai) + C 2 (//) + 6C 3 (/i) + 60C 5 (//), (D.5) 

7 2 32 

= ~2 C M ~ - 38C M - yCe^), (D.6) 

i 9 32 

w 9 = --C 3 (/i) - -cm - zcm - -CM, (D.7) 
where C;(/x) = cf V) + a s (n) / (li^C? \n) + • • • and the 7 J 9 0) 's can be found in Ref. [5]. 
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